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Abstract 

A mathematical formalism is presented which allows studies in multistep processes where 
one particular group of objects is generated and in turn generates a new group with similar 
properties with some rates. Such multistep processes are observed for example in epidemiology. 

For this purpose an adaptation is made of the Bateman equations describing a system of 
nuclear levels interconnected by transitions. Closed mathematical expressions allow an easy 
calculation/prediction of the properties of such multistep processes. 

The need of adequate mathematical description/modeling of epidemics evolution appeals strong¬ 
ly for developments in this field, it is enough mentioning the pandemic COVER-19 in the beginning 
of 2020 (see e.g. the report[l] of the WHO.) The author of the present work is an amateur in epi¬ 
demiology, but a brief consideration of recent papers in that field employ a variety of mathematical 
apparatus to analyze the statistics of the available data in order to provide explanations for it and 
predict possible future scenarios (cf. e.g. Refs. [2, 3, 4, 5] and references therein). While an analysis 
of the initial part of the experimental curve describing the number of observed cases as function 
of time is relatively reliable for short range forecasts, also when the forecast is very close to its 
maximum (peak), in general it is desirable to dispose with description of the epidemics in the full 
time range till its disappearing (falling in the background). More involved epidemics models [6] 
involving also solving of non-linear differential equations (e.g. SIR) are given in Refs. [7, 8]. 

One different possibility for epidemics description is provided by adapting the Bateman equa¬ 
tions known since the beginning of the 20th century (see e.g.[9], also [10] for more details), which 
are linear, for the epidemiology case. The Bateman equations describe the time evolution of the 
population of a system of nuclear levels which are interconnected by transitions (e.g. emission of 
gamma-rays in Nuclear physics, also photons in similar systems in Atomic physics). Namely, the 
system of N nuclear levels is governed by N differential equations (i=l...N): 

hi(t) = -X+ Y X h b hi n h (t) (1) 

h 

where n*(t) is the evolution of the population of the level of interest i, the levels h (h > i) are 
direct feeders of that level, Aare the decay constants and bf n are the branching ratios. The 
general solution (if Xk ^ X j for i ^ k) including all possible feeding levels can be written in a 
straightforward way (c.f. e.g. ref. [9]) as 

n,{t) = Y, C ik e~ Xkt = (m( 0) - Y Cik)e~ Xit + Y C ik e~ Xkt (2) 

k>i k>i k>i 

where n,(0) is the initial population (at t = 0) of the level i, the numbering of the levels starts 
from it and the coefficients are given by the recursive equations (k > i) 

Cik = t— (bkiHk( 0) - Ckmbki + Y'' bu^Cik] (3) 

* k V m>k l>i k / 

On the other hand, provided the functions nh(t) of the direct feeders h are known the solution 
of the set of equations (1) for every particular level i could be written as: 



In the latter formula, no details about the history of the population of the feeding levels h are 
present. However, this information is contained in the functions rih(t) and that is the way the 
function n,(f) ’’feels” the earlier stages of the cascade history. 

We use the above approach to describe in very crude way a process which has the features of 
epidemics. Only two decay rates Xr (rate of recovery) and Ac (rate of contamination) are employed 
which are assumed to be constant in time. Of course, the latter may be not true in the history of 
each case of disease when going from the begin of the illness, the peak (crisis) and the recovery 
(see discussion below). Let us consider the evolution of the different generations of contaminated 
people with time. At time t=0, there are ni(0) cases which recover with a constant rate Xr. 
In practice, such assumption may be very crude since: (i) the time necessary for the recovery 
process depends on the health parameters of the patient in the more general sense, (ii) the rate 
of recovery may be also not a constant i.e. the illness goes to different stages during the recovery, 
characterized by different health parameters (iii) application of treatments with medicaments with 
different efficiencies. Nevertheless, an average of A^’s can be used in calculations as done below. 
Also, some upper limit A r fixes the time limits of the full process of acquiring immunity for the 
whole population. The same considerations apply to the rate of contamination Ac which, however, 
may be very well influenced by measures as confinement of ill people, reduction of social contacts 
etc. With these remarks the time evolution of n\(t) is governed by the well known radioactivity 
decay law of Rutherford-Soddy: 

^ = -W) (5) 

with the solution n\ ( t ) = n\ (0)e~ Aiit . Here, the mutual new contamination of already contaminated 
people is not taken into account, of course. 

The population of the second generation with time 712 (f) is governed by the equation 

dU ^ = -A R n 2 (t) + X c ni{0)e~ XRt (6) 

It is easy to show that its solution, with the boundary condition n2(0)=0 is given by: 

n 2 (t) = A c ni(0)fe _ABt (7) 


The third generation population 77 . 3 (f) is governed by the equation 

dn 3 (t ) 


dt 


= -X R n 3 (t) + X c n 2 (t) 


Again it is easy to show that its solution is given by: 

Ac 2 ni(0)f 2 e“ AKt 
Mt) = -o-> 


( 8 ) 

(9) 


for the fourth generation one obtains 714 (f) = Ac 3 7 ii( 0 )f 3 e Ah 4 /(2.3) and so on. The final result 
for the n th generation obtained by induction reads: 


n n {t) 


Ac"" 1 m(0)f"~ 1 e- ARt 

(n- 1 )! 


( 10 ) 


The summing up the number of contaminated people from different generations one obtains for 
their total number as function of time: 


N(t) 


E 

k—1,2,..NMax 


(^) k ~ 1 ni{0)x k ~ 1 e~ x 


(k - 1 )! 


( 11 ) 


where Nmux is the last generation and x = A Rt. This expression represents a sum of Poisson 
distributions weighted by the factors (y^) fc_1 - 



Table 1: Calculations for contamination rate A c = 20 contaminations per day caused by an ill 
person and different mean recovery times tr (Xr= \/tr). The initial population ni(0) = 1. In 
the second column, the the contaminated generation number n is shown followed by the position 
of the peak maximum in time tMax• For the other quantities see text. The table is left-aligned 
and no abbreviated formats for the large numbers (as e.g. a*10 b ) are used in order to get better 
optical effect on the observed increase with n and tr. 


tr 

n 

tMax [days] 

(n — l)" _2 e _ l n_1 l/(n — 2)! 

(Xc/Xr)^ 

^Max 

n n 

7 days 

2 

7.0 

0.368 

141 

52 


3 

14.1 

0.135 

19837 

2684 


4 

21.1 

0.075 

2793990 

208656 


5 

28.2 

0.049 

393519814 

19220178 


6 

35.2 

0.035 

55425326023 

1945067235 

14 days 

2 

14.1 

0.368 

281 

103 


3 

28.2 

0.135 

79349 

10738 


4 

42.3 

0.075 

22351925 

1669255 


5 

56.3 

0.049 

6296317036 

307522851 


6 

70.4 

0.035 

1773610432745 

62242151525 

21 days 

2 

20.8 

0.368 

417 

153 


3 

41.7 

0.135 

173611 

23496 


4 

62.5 

0.075 

72337962 

5402242 


5 

83.3 

0.049 

30140817901 

1472128897 


6 

104.2 

0.035 

12558674125514 

440727503315 


It can be shown that the function n n (t) has one maximum which occurs at 

tmax = {n-l)/\ R = (n-l)T R (12) 

where with tr we denote the reciprocal of A r and which is equal to the mean time needed for 
recovery (analog of the mean lifetime of an excited nuclear state). Thus, at tr = 14 days the 4 th 
generation of contaminated people will have its maximum (peak) 42 days after the begin of the 
contamination, the 5 th one - after 70 days and so on within our simple model. The height of the 
peak maximum at t max is given by 

n n (t) = (^)"- 1 n 1 (0) (n ~ 1) "~ 2 e- (,| - 1) . (13) 

A r (n-2)1 

and a more careful consideration (see also Table 1) reveals that its value depends extremely strong 
on the ratio {Xc/Xr) at the power of (n— 1). For example, if that ratio is 1000, for the 4 th generation 
the result will lead to 10 9 . The factor (n — 1 ) n ~ 2 e~ < - n ~ 1 ^ > /(n — 2)! will reduce that number only 
by about 10%. Another interesting feature of the formalism is that in the composition of N(t), 
the total number of contaminated people, the last generation of contaminated riN max (t) dominates 
also strongly the other terms of the sum (with a factor of 100 and more, Table 1) at any time. 
Then, the peak position of n,N m ax(t) dominates the peak position of the development in time of the 
whole epidemics, indeed. Some more details, e.g. on the time range of the epidemics, are presented 
in Fig. 1 

This paper does not intend to investigate and provide predictions for the development of con¬ 
crete epidemics and pandemics. Such tasks are beyond its scope and the competence of the author. 
But our results may be used by specialists in epidemiology when the remarks made above about 
the rates A r, Xc are taken into account. Moreover, one has to consider important factors as 
specificities in every country coexisting with globalization effects, the existence of different social 
groups and different health systems in each country etc. From that point of view, we think that 
the mathematical formalism proposed here can help the specialists in epidemiology to solve prob¬ 
lems and if necessary, they can introduce modifications in the approach to take into account more 
complex situations. 





Figure 1: A graphical representation of n n (t) functions of third an fourth generation combined in 
pairs with common tr (7 or 14 days) and a common Ac=20 for all four cases. Note the different 
factors used to scale apart that of tr= 14, n=4. The time extent of the epidemics for tr=7 and 
tr= 14 for the cases displayed is about 100 and 200 days, respectively. See also text. 
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